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Novel Properties of Bound States of Klein-Gordon 
Equation in Gravitational Field of Massive Point 

p.p. Fizievt T.L. BojadjievJ D.A. Georgieva ^ 



Abstract 



We consider for the first time the solutions of Klein-Gordon equation in grav- 
^sq ' itational field of massive point source in GR. We examine numerically the basic 

bounded quantum state and the next few states in the discrete spectrum for dif- 
^^ I ferent values of the orbital momentum. A novel feature of the solutions under 

\^ • consideration is the essential dependence if their physical properties on the gravita- 

cn . tional mass defect of the point source. Such mass defect was not introduced up to 

recently. Its variation yields a repulsion or an attraction of the quantum levels up 
Q ■ to their quasi-crossing. 

^' 

o 
^' 1 Introduction 

^1 The correct solutions of the Einstein equations for the gravitational field of point particle 

with bare mechanical mass Mq were found recently H^. It turns out that these fundamental 
solutions are described by mathematical distributions and own the necessary jumps in the 
first derivatives of the metric, needed to satisfy the Einstein equations at the place of the 
point source. The new solutions form a two-parameters family of metrics on singular 
manifolds Wl^^'^^{g^^}. They are defined by the bare mass Mq and by the Keplerian mass 
M, or equivalently, by the Keplerian mass M and the ratio of the masses: g = M/Mq G 
(0, 1). This ratio describes the gravitational mass defect of the point particle. 

The mathematical and the physical properties of the new solutions are essentially 
different in comparison with the well known other spherically symmetric static solutions 
to the Einstein equations with different type of singularities at the center of the sym- 
metry, which is surrounded by an empty space. The previously known solutions were 
often considered as a solutions for describing of single point particle in GR. In the most 
of the known cases this turns to be incorrect ^J. Therefore the new solutions call for 
reconsideration of many of commonly accepted notions and their physical interpretation. 

Here we are starting the study of the GR wave equations in the gravitational field of 
massive point. In the present article we shall consider as an example only the simplest 
GR wave equation - the Klein-Gordon one. We shall present a numerical study of this 
well known equation which describes (anti)particles of zero spin jHj. 
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The reader can find the basic resuhs of the analysis of this equation in the vacuum 
Schwarzschild metric, when described in the commonly used Hilbert gauge or in other 
gauges, for example, in the recent review article |3], in j3], and in the large amount of the 
references therein. A discussion of different radial gauges and corresponding references 
one can find in pp. 

2 The Formulation of the Problem 

Making use of the natural radial variable of the problem: r G (0, cxo) P], one can represent 
the new regular solutions for the gravitational field outside of the massive point source in 
the form: 



Here 



df ^' 



Na{ry 



-p{rf{de^+siYi^ed(f). (2.1) 



is the modified Newton potential, the coefficient N^[r) = {2ip^)~ (e^'^c — 1) , the Hilbert 
luminosity variable is 

_ 2G^M _ r + G^M/H^) 

and Gj^ is the Newton gravitational constant. 

In Hilbert gauge, outside of the source, the solution has a standard form: 

gtt{p) = l-pG/p, Qppip) = -l/gtt{p), (2.4) 

where pc = 2G^M is the Schwarzschild radius. We have to stress that the presence of 
the matter source forces us to consider this form of the solution only on the physical 
interval of the luminosity variable p ^ {po, oo), where 

po = 2G^M/(1 - g^) > pG. (2.5) 

One can find in [T^ the physical and the mathematical justification of this cutting pro- 
cedure of the admissible interval of luminosity variable p for the discovered by Hilbert, 
Droste and Weyl |2J form of the solution ()2.4|) . Here we shall stress that such unusual 
cutting is a ultimate consequence of the existence of a matter source of the gravitational 
field. After all, it permit us to satisfy the physical requirement not to allow an infinite 
luminosity of the point sources in correspondence to the real observations. In the articles 
PP it was shown that in GR we have a natural cutting parameter for the classical diver- 
gences. It originates from the drastic change of the space-time geometry due to the infinite 
mass-energy density of the point particle. It turns out that in presence of matter source 
all phenomena, related with the event horizon, together with the very horizon, belong to 
the non-physical domain of the variables. Therefore, in accord to Dirac's intuition jH], all 
such notions "should not be taken into account in any physical theory" of matter. 



The luminosity variable p is an exclusive one, because the results, expressed by making 
use of it, are locally invariant under radial gauge transformations. From computational 
and from physical point of view a more preferable variable is: g = gu = I ~ Pg/p- Then 



ds"^ = g dt^ — pQ 



dg^ de^ + sin^edcP^ 



,g{l-gr (l-gy 



(2.6) 



In these variables the 4D D'Alembert operator in the curved space-time M'^^'^-'j^f^,^} 



fii^j ) 



created by a massive point, acquires the form 

□ = g-'d^ - Pg' ((1 - 9)% iad,) + (1 - g^Ae^) (2.7) 

where Ag^ = sm~^9dg (sinOde) + sin~^6'9^. 

The Klein-Gordon equation for test particles of mass m in space-time Wl^^'^^{g^^} has 
a standard form: 

n$ + m2$ = 0. (2.8) 

As a result of the spherical symmetry the pseudo-Reimannian manifold M.^^'^'{g^^}, 
created by the point particle, has a group of motion SO (3). Nontrivial are the quantities 
and equations only on the quotient space M^^'^^ = Wl^^'^^ /S0{3), with natural global 
coordinates t and g. The reduction of Eq. ()2.8|1 on the quotient space M*-^'^^ can be reached 
by the substitution ^{t,g,9,(f)) = '^i(t,g)Yi^i^{9,(f)), where Yi^i^{9,(f)) are the standard 
spherical functions: Ag^Yi^i^{9,(j)) = -l{l + l)Yij^{9,(l)), I = 0,1,2,...; 4 = -/,..., 0, ..., Z. 

For the time- dependent radial wave function '^i(t,g) on the manifold M*^^'^) we have 
the partial differential equation: 

g-'d^^i - p^'{{l - g)% {gdg^i) +[/(/ + 1)(1 - gf + m"] ^i] = 0. (2.9) 

The space-time M'-^''^^{(7^,y} of the problem at hand is invariant under translations 
with respect to the time variable t, as well. This symmetry of M.^^'^'{g^^} is related to the 
translation group T(l): t —>■ t + const and allows the total reduction of the problem to one 
dimensional one using the anzatz '$i{t,g) = c~'^^^Ri{g). For the time-independent radial 
wave function Ri{g) on the quotient space M^^^ = M(^'^Y^(1) 'we obtain an ordinary 
differential equation: 



d^Ri 1 dRi 

H + 

dg"^ g dg 



e" p" 1(1 + 1] 



,g\l-gy g{l-gy g{l 



Ri = 0. (2.10) 



Here e = pc E and p = Qg^^ are, correspondingly, the dimensionless total energy and 
mass of the spinless test particles. We are using the units c = h = 1. The dimensionless 
orbital momentum in our units is / = L/mpc, where L is the dimensionfull one. 

The motion of the Klein-Gordon particles with a fixed orbital momentum / can be 
treated as a relativistic motion in a dimensionless potential 

vi{g) = (1 - (?) [/(/ + 1)^(1 -g)-p'], ge [go, 1], (2.11) 

which is shown (for p = 1 and several different values / > 0) in Fig. |21 The point g = 
corresponds to the event horizon and is placed in the nonphysical domain g G [0,go), 
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Figure 1: The potential f;(5') for / = 0, 1, 2, 3, 4 and /i = 1. The dashed hne corresponds 
to the value /^"* = |a/1 + 12/1^ — | for which the maximum of the function vi{g) arise. 



where f^o = f?^ > 0. The point g = ^ corresponds to the physical infinity (with respect 
to the variables r, or p). There the space-time M'^^'^^j^f^i,} is flat. Writing down the Eq. 
flTTmi in the form 

d 



(1-#U-J Ri^\\-vi{g)\Ri = ^ 

we see that the points (7 = g = 1 are singular points of this equation. The second one 
is a non-regular singular point. Here A = e^ — /x^. 
The substitutions: 

' g \ 1 



M = In 



i~gj {i-gY 

Ri{g) = {l-g)Pi{g), 



(2.12) 



introduce "the tortoise coordinate" u and the corresponding radial wave function Pi{u) 
and lead to the following standard (Schrodinger like) form of the Eq. ()2.10|) : 






+ [\-wi{u)\Pi = Q^ 



(2.13) 



The bounded quantum states can be found as a solutions of ()2.13|) under the boundary 
conditions 



together with the usual L2-normalization 



(2.14) 



oo 

/ 



uo 



Pfiu)du- 1 = 0. 



(2.15) 



The quantity 

wi{g) = vi{g) + g{l - gf = {I - g) [l{l + 1)^7(1 ~ g) + g{l - gf - ^l^] (2.16) 



is the extended dimensionless potential of the problem and depends implicitly on the 
variable u via the solution ()2.12p of the Cauchy problem 



du 



5'(l-5')^ 9^ = 0^ 



(2.17) 



An essential novel feature of our approach to the problem is the correct fixing of the 
physical domain of the variables: the potentials vi and wi must be considered only on 
the interval [uq.co), where Uq = Uq{q^). Thus the nonphysical infinitely deep potential 
well at p ^ 0, and the event horizon at p = pg{^ug = —oo), are excluded from our 
consideration. 
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Case: Q^ = 0.\,n = 
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Figure 2: The radial wave function Poi{u, q^) for different values of the quantum number 



/ = 2, 3, 4, 5, 6, 7,: 



As a completely new phenomenon, a discrete spectrum comes into being in the problem 
at hand. We denote hy Sni = yKii + fJ^ the discrete values of the energy e of the 
bounded states. The principle quantum number n = 0, 1, ... determines the number of the 
zeros of the radial wave function Pniiu) in the interval u G (mq, oo) and of the function 
Rni{g) = (1 - g)Pni{u{g)) - in the interval g G (5-0, 1)- 



3 Basic Numerical Results 

For the numerical study of the solutions of the Sturm-Liouville problem ()2.13p - ()2.15|) 
on some finite interval [mq, Wqo] we use the continuous analog of Newton method [Xl- The 
corresponding linear boundary problems at each iteration are solved by means of spline- 
collocation scheme of fourth order of approximation. The convergence of the calculated 
eigenfunctions and corresponding eigenvalues A„; for large enough values of Mqo is proved 
on a set of embedded intervals. As a result of numerical experiments we have obtained a 
number of new properties of the bounded quantum states of the problem at hand, a part 
of which were unexpected to some extend. 

The basic results for the radial wave function Pni{u) for different values of the orbital 
quantum number / an n = are shown in Fig . |21 
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Figure 3: The radial wave function Po2{u; g^) for different values of the two variables u 
and the squared mass ratio q^. 



Case: / = 2, « = 




Figure 4: The radial functions Po2(w; Q^) for different values of the gravitational mass 
defect of the source. 



In the Fig. El using as an example the wave function Po2(w; ^^), we show the typical 
dependence of the wave functions of the bounded states both on the variables u and on 
the squared mass ratio g^. 

The radial functions Po2('w; Q^) for different values of the gravitational mass defect g 
are shown in Fig. |3] As seen, depending on the values of the squared mass ratio g^ of the 
point source, the wave function of the test particle in the point particle's gravitational 
field is located in the inner potential well (for small g^ - the curves 1 and 2 in Fig. |3]), 
or in the exterior potential well (for big g^ - the curve 4 in Fig. H}. In the narrow 
intermediate domain of values of g^ a transition regime take place - the curve 3 in Fig. EJ 
Actually, this transition depends continuously on the variable g^, but it seems that it is 
a jump-like, because the transition happens in a quite smaller scales with comparison to 
the other much more smooth changes of the corresponding quantities. (See, for example, 
the dependence of the eigenvalues on the variable g^). These phenomena are described in 



a more transparent way in the 3D Fig. IHl 

In our approach to the problem the basic new physical phenomena are related with the 
gravitational mass defect of the point source. Such mass defect was not considered and 
studied until now, because for the standard Hilbert form of the Schwarzschild solution 
"the bare rest-mass density is never even introduced" [Hj correctly. 

The dependence of the first four eigenvalues \o 1,11,21,31 on the squared mass ratio g"^ 
is shown in Fig. El for / = 2. A typical steep-like dependence of the discrete eigenvalues 
^niiQ"^) on the variable g"^ is seen. For each value of the principle quantum number n the 
number of the jumps of the energy levels Sniig^) = \/^ni{g^) + /W^ depends on the number 
of the maxima of the corresponding wave function which are moved from the inner well 
to the exterior one during the transition regime which develops with the increase of the 
values of the mass ratio g^. 
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Figure 5: The dependence of the first four eigenvalues \oi,ii,2i,3i on the mass ratio g"^. 



It's clear that such nontrivial behavior of the energy levels eni{g'^) of the Klein-Gordon 
equation in a gravitational field of a point particle is due to the existence of two finite 
potential wells: 1) An inner one which has a size of order of the Schwarzschild radius 
and in general case of macroscopic orbital momenta is very deep; and 2) An exterior 
one, which is extremely shallow in comparison with the inner well. The normal world 



with almost Newton gravity "lives" in the exterior well. The both wells are separated 
typically by a huge potential barrier which looks like a centrifugal barrier from outside, 
and as a suspensory barrier of the type of a spherical potential wall from inside. Our 
calculations describe the quantum penetration under this barrier. The quantum result 
depends strongly on the mass defect of the point source. 

The obtained in this article behavior of a quantum test particles in the gravitational 
field of point source of gravity seems to us to be much more physical then the one in 
the wide spread models of black holes. Clearly, in contrast to such space-time holes with 
nonphysical infinitely deep well in them, in our case the finite inner well plays the role 
of a trap for the test particles. It is not excluded that this way one may be able to 
construct a model of very compact matter objects with an arbitrary large mass and a 
size of order of Schwarzschild radius. It's possible that such type of objects may explain 
the observed in the Nature compact dark objects. They may be the final product of the 
stelar evolution, instead of the quite formal and sophisticated constructions like black 
holes which correspond to empty space solutions of Einstein equations. Of course, at 
present these possibilities are an open problem which deserves further careful study. 
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Figure 6: The attraction and the repulsion of the quantum energy levels for different 
values of the gravitational mass defect of the point source. 



It is easy to observe one more astonishing phenomena in the discrete spectrum of a 
test quantum particle in gravitational field of point source. It is related, too, with the 
mass defect of the source. If one puts in the same figure the functional dependence of 



all squared discrete energy levels 5„; on the squared mass ratio q 



9tt 



see m 



E), 



one can observe a repulsion and an attraction (up to a quasi-crossing) of these levels, 
as shown in Fig. Such type of behavior of quantum discrete levels is well known in 
the laser physics, in the neutrino oscillations and in other branches of quantum physics. 
To the best of our knowledge we are observing such behavior of quantum levels in the 
fundamental gravitational physics for the first time. 



4 Conclusions 

In the present-days traditional approach to the motion of test particles in the 
Schwarzschild gravitational field, inside the event horizon there exist a nonphysical in- 
finitely deep potential well. Everything which somehow can be placed inside this well will 
fall unavoidable to its center p = for a finite time. It is well known that this center is 
a physically inadmissible geometrical singularity. To hide such undesired singularity one 
must use a unprovable mathematical hypothesis like the cosmic censorship one. Today 
we have enough examples which show that such hypotheses can not be correct at least in 
its original formulation, see, for example [9^ and the references therein. 

In the article pi the principal role of the matter point source of gravity was stressed. 
In GR the matter point source presents a natural cutting factor for the physical values of 
the luminosity variable p G [po,oo), where po > Pg fj^-Sj) . This is because the infinite mass 
density of the matter point changes drastically the geometry of the space-time around it. 
In full accord with Dirac's suggestion |^ this cutting places all nonphysical phenomena, 
together with the very event horizon in the nonphysical domain of the variables and yields 
a large number of new interesting phenomena. 

In particular, as shown in the present article, the described cutting yields an interesting 
novel discrete spectra for Klein-Gordon test particles in the gravitational field of massive 
point. Mathematically this is caused by the universal change of the boundary conditions 
for the corresponding wave equation, due to the presence of the matter source of gravity 
and depends on its mass defect. 

A real discrete spectrum does not exist in the case of black hole solution, just because 
of the presence of a hole in the space-time, which absorbs everything around it. Our 
consideration gives a unique possibility for a direct experimental test of the existence of 
space-time holes. If one will find a discrete spectra of corresponding stationary waves, 
propagating around some compact dark object, one will have indisputable evidence that 
there is no space-time hole inside such object. If, in contrast, one will find only a decaying 
quasi- normal modes with complex eigenvalues (see IS] and the references herein), this w 
ill indicate that we are observing indeed a black hole. 

We have obtained similar results for electromagnetic and gravitational waves in the 
gravitational field of point particle ^U]- For them the properly modified quasi- normal 
modes turn to depend on gravitational mass defect, too. The further study of these phe- 
nomena is an important physical problem and will give us more practical physical criteria 
for an experimental distinguishing of the two complete different hypothetical objects: 1) 
the space-time black holes and 2) the possible new very compact dark objects made of real 
matter. Present days astrophysical observations are still not able do make a difference 
between them. The only real observational fact is that we are seen very compact and very 
massive objects, which show up only due to their strong gravitational field. An actual 
theoretical problem is to find a convincing model for description of these already observed 
compact dark objects and the criteria for the experimental justification of such model. 
We hope that the present article is a new step in this direction. 
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